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Abstract 

The recent experimental data on D + — D° and D* + — D*° mass differences are used as inputs 
in the QCD sum rules to obtain new estimates on the mass difference of light quarks and on 
the difference of their condensates: vn^ — m u = 3±1 MeV , (dd) — (uu) = — (2.5±l)-10~ 3 {uu) 
(at a standard normalization point, fi = 0.5 GeV). 
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The QCD sum rules invented more than a decade ago is now well known to be a very 
useful tool to study properties of hadrons at intermediate energies and to get information 
on the basic parameters of QCD, such as quark masses and non-perturbative condensates. 
One of the problems addressed already in the pioneering papers[|IJ was the relation between 
the isotopic symmetry violation on the level of hadrons and the difference between u- and 
d-quark masses and condensates. It was shown that the observed p — uj mixing implies 
( m.<2 — m u ) j (m u + m<i) ~ 0.3 and 7 ~ —1.5 • 10 -2 , where 

7 = (dd)/(uu) - 1 (1) 

ruling out a solution with m u = 0, 7^ 0. The quark mass difference was estimated 
before the advent of QCD sum rules in refs. |2|, |3|] with the result — m u « 3 MeV and 
fri'd J r m u ~ 11 MeV . The difference of condensates was later estimated within the QCD sum 
rule framework in a number of papers^] with the result 7 = — (3 4- 10) • 10~ 3 . The difference 
of condensates was also obtained in the framework of chiral perturbation theory 0, which 
gives 7 ~ —8 • 10~ 2 , provided — m u = 3 MeV, 1 — (ss) / (uu) ~ 0.2 and m s m 150 MeV. 
Recently isospin violation in QCD sum rules for the nucleon, S and 5 was considered || and 
the following results were obtained: — m u = 3 ± 1 MeV, 7 = — (2 ± 1) ■ 10~ 3 . Thus, 
while most predictions for — m u agree and are grouped around 3 MeV, predictions for 
(dd) — (uu) are more diverse and range within an order of magnitude. Moreover, arguments 
were given in ref.0 that m u may be equal to zero due to instanton contributions to the 
renormalization of the quark mass. Thus, additional independent estimates of differences 
between u— and d— quark masses and condensates are certainly welcome. 

In this paper we will consider QCD sum rules for isospin mass splittings of D* and D 
mesons and make use of the recently reported f| new results on these splittings, 

m D «+-m D *o = 3.32 ±0.08 ±0.05 MeV 
m D + - m D o = 4.80 ± 0.10 ± 0.06 MeV (2) 

to obtain such estimates. The sum rules are similar to those which were used in ref.0 
to succesfully predict the mass splittings m^* — mo* = HO ± 20Mel^ and mo s — trip = 
120 ± 20 MeV. 

Let us start with the correlator of two pseudoscalar currents with quantum numbers of 
D, J5 = c-faq, where q is either u, or d, at Euclidean momentum — q 2 > 1 GeV 2 , 

Ci = tJ rf 4 xe^(0|T{j 5 (x), j 5 + (0)|0) (3) 
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and consider its variation 5C q as the light quark mass rises from zero to its actual value 
m q . To estimate C q it is possible to take into account only the unit operator and the quark 
condensate (Fig.l) in the operator product expansion, since the contribution of operators of 
higher dimension to the heavy-light correlators is negligible|TD|. Using the expansion of the 
quark condensate in the quark mass 

(Qa^Hm = (~Sap(qq) + ^m q x a p(qq)^ 5 ab (4) 

where (qq) itself also depends on m q , it is easy to see that 5C q is a function of m q and 
(qq) — (qq)o where the subsript denotes the chiral limit. On the other hand, saturating 
the correlators in the standard manner by the corresponding lowest mass resonanses D + and 
D°, subtracting the continuum from the contribution of the unit operator and applying the 
Borel transformation [||], (s + Q 2 )^ 1 — > M~ 2 exp(— s/M 2 ), we arrive at the following sum rule 
for 5C d - 5C U 



(3d+-(3d°^2 , x M A e</ M ' 2 



M + {m D + — m D Q 



Pd^D "' ' V " ±J ' ■""~'hadr 2(3' D m D 

{(m d - mj[^(e- - e"» - x(E 1 (x) - E 1 (y)))L^ - ^e~*(l + ^)] 
^Juu) ^^ _ {sd+ _ Sdo) _±_ {sd _ mlfe -y } (5 ) 

Here E\(x) = f£° dtt~ 1 e~ t , m c is the charmed quark mass, sd is the continuum threshold, x = 
m 2 /M 2 , y = sp/M 2 and the residue of D meson into the current j'5 is f3 D = /B m i)/ m c) where 
fo is the semileptonic decay constant defined by (_D| 07^755 10) = —ifoV^- The dependence 
of light quark masses and condensates on the normalization point in the operator product 
expansion is given by powers of L = ln(M/A)/ ln(///A) where A = 150 MeV and we take 
= 0.5 GeV which corresponds to (qq) = -(0.24 GeVf. Here (3 D + = f3 D + 513%, f3 D o = 
Pd + 5(3 D and s D + = sn + 8s D , s D o = sn + 5s D where (3d and sd are the residue and the 
continuum thresholds in the chiral limit and 5j3 D ,d and 5s u D ,d are the deviations of residues 
and of continuum thresholds from their values in the chiral limit. 

Similar sum rules can be written for the correlator of two vector currents, j M = c^^q, 

C% = if d*xe iq *(0\T{j,(x), j+(0)|0> (6) 



In the case of the tensor structure q^q u for which the sum rule is known to work better [|Tl~1 
we obtain 



2 



,Pd*+-Pd*° (m D *+ - m D *o) hadr 2 

)M + (m D .+ - m D *o) hadr 



f3 D *m D * m 2 D 

m 2 //l/f 2 

rriD*e 



5 [(m d - m u )(gg)e + (s D ,+ - s D *o)— —^-{l - + -^)] (7) 



It is important to emphasize that since we do not take into account perturbative two- 
loop diagrams with one photon exchange in the "theoretical" part of the sum rules, only 
the hadronic parts of the isospin splittings, (Am D ) hadr and (Amo.)u n enter eqs.flU) and 
(0). To obtain them we use a quark model estimate [0] for the photon cloud part of the 



mass difference m D + —m D a = 1.7 ±0.5 MeV and also take into account the electromagnetic 
hyperfine splitting 

Sm = j2 1 ommi (8 ) 

3m c m g 

where S is the spin of the meson. The quark masses here are the constituent ones, m c ~ 
1.7 GeV, m q ~ 0.3 GeV. The quark- ant iquark wave function at the origin may be estimated 
using the relation /|, = 12|\E , (0)| 2 /m £) . Using the estimate flCfl fr> = 170MeV, we get that the 



hyperfine electromagnetic interaction contributes ~ 0.5 MeV to D + — D° and ~ —0.17 MeV 
to D* + — D*° mass splittings. Using the full experimental mass differences from eq.(|2]) we 
thus have the mass splittings to be used in the sum rules in eqs.(^) and (|7|) 



(%+-%•«)/.* = 1-8 ±0.5 MeV 
{m D+ - m D o) hadr = 2.6 ±0.5 MeV (9) 

From the sum rules in eqs.(||]) and (0) it follows that in the working region for the Borel 
parameter M 2 , the left hand sides should be close to linear functions in M 2 whose extrap- 
olations to M 2 = give the corresponding hadronic mass splittings while the slopes give 
information on the residue differences, Pd+—Pd° and /?£>. + — Pd*°- To numerically analize the 
sum rules we take mjj = 1.87 GeV, m* D = 2.01 GeV, the standard value of the quark conden- 
sate, (qq) = — (0.24) 3 GeV 3 , and m c = 1.35 GeV. For the residues and continuum thresholds 
we take the estimates obtained with sum rules in refs. |T0|, 11| , /3f, ~ /3|>, ~ 0.2 GeV 4 and 
s D « s D * « 6 GeV 4 . 

Thus, m d — m u , (dd) — (uu), the Borel parameter M 2 and the differences of continuum 
thresholds, sd+— sd° and Sd* +— Sd*o are the fitting parameters of the sum rules in eqs.(^|) and 
(0). The working regions in M 2 determined by the requirement of controllable contributions 
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from non-perturbative power corrections and continuum in the chiral limit for u- and d- 
quarks are[|l(], [TIJ] lGeV 2 < M 2 < 2GeV 2 for the pseudoscalar channel and 1.5 GeV 2 < M 2 

< 2.5 GeV 2 for the vector channel. We will let the difference of thresholds vary around an 
estimate Sd+ — s^o ~ Se>*+ — So*o ~ (m d — ?ti m ) v /Jd which implies that the threshold and 
the quark mass differences are of the same sign. 

In Fig. 2 we show the result of the sum rule calculation of the D* mass difference by 
plotting f v {M 2 ) - M 2 df v (M 2 )/dM 2 , where f v {M 2 ) is the r.h.s. of eq.(0). From eq.® 
one sees that only the quark mass difference enters this sum rule and one might hope to 
fix m d — m u from it. The term proportional to S £)* + — S £)* has the same sign as the term 
proportional to m d — m u and in the whole working interval in M 2 contributes less than 50% 
to f{M 2 ) - M 2 df(M 2 ) jdM 2 at m d - m u = 3 MeV. The values m d - m u > 4 MeV cannot be 
made compatible with {m D *+ — m D *o) hadr from eq.(|9|) at any Sd*+ — Sd*° > 0. At m d — m u 

< 4 MeV the compatibility can be achieved by varying sd*+ — sn*o. However, at m d — m u 

< 2 MeV the continuum term is dominant in the whole working region in M 2 and the sum 
rule is not reliable. 

As for the D mass difference, the corresponding sum rule in eq.(^) is contributed both 
by m d — m u and (dd) — (uu). In Fig. 3 we show the numerical results for m d — m u = 3 MeV 
and two different values of 7 = (dd) / (uu) — 1. One can see that the value 7 = —2.5 ■ 10~ 3 
is consistent with the (m£>+ — m]jo) hadr from eq.(||), while a higher value 7 = — 6 • 1CT 3 is 
definitely excluded. 

In summary, we conclude from our numerical analysis of the sum rules for isospin mass 
splittings in D and D* mesons that m d — m u = 3 ± 1 MeV and 7 = —(2.5 ± 1) ■ 10~ 3 . The 
result for m d — m u is consistent with the earlier estimates and together with the relation 
(m u + m d )(uu) = —F 2 ml excludes the option m u = advocated in ref. 0. Our estimate for 
the condensate difference parameter 7 supports the value obtained in ref. from the isospin 
splittings in baryons and differs from results of earlier papers^, [JQ. 

We are grateful to H. Leutwyler for useful discussions. One of the authors (V.L.E.) 
acknowledges the warm hospitality extended to him at the Universities of Minnesota and 
Bern. This work was supported in part by the US Department of Energy under contract 
DE-FG02-87ER40328 and by Schweizerischer Nationalfonds. 



1 For a discussion of possible reasons for this disagreement, see ref. [ 
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Figure Captions: 



• Fig. 1: The diagrams taken into account in the sum rules. 

• Fig. 2: Sum rule for D* mass splitting: f v (M 2 ) - M 2 df v {M 2 ) / dM 2 versus M 2 , where 
fv(M 2 ) is the r.h.s. of eq.(0) for m d -m u = 1 MeV (a), 3 MeV (b) and 5 MeV (c). In 
each of the three cases the lower and the upper curves correspond to so*+ — sd*» = 
and 0.005 GeV 2 , respectively. The dotted horizontal lines are the boundaries set by 
eq.@. 

• Fig. 3: Sum rule for D mass splitting: f P (M 2 ) — M 2 df P (M 2 ) / dM 2 versus M 2 , where 
f P {M 2 ) is the r.h.s. of eq.flD for m d - m u = 3 MeV and a) 7 = -2.5 ■ 10~ 3 ; b) 
7 = — 6 ■ 10~ 3 . The numbers at the curves correspond to the value of sd+ — s^o (in 
GeV 2 ). The dotted horizontal lines are the boundaries set by eq.(|]). 
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